Exact Solution of a Drop-push Model for Percolation 
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Motivated by a computer science algorithm known as 'linear probing with hashing' we study 
a new type of percolation model whose basic features include a sequential 'dropping' of particles 
on a substrate followed by their transport via a 'pushing' mechanism. Our exact solution in one 
dimension shows that, unlike the ordinary random percolation model, the drop-push model has 
nontrivial spatial correlations generated by the dynamics itself. The critical exponents in the drop- 
push model are also different from that of the ordinary percolation. The relevance of our results to 
computer science is pointed out. 



PACS numbers: 64.60.A, 02.50.-r, 05.40.-a, 89.20.-a 

The ordinary site or bond percolation and its vari- 
ous generalizations are amongst the most well studied 
problems in statistical physics ||,|2j . Motivated by a well 
known computer science algorithm known as 'linear prob- 
ing with hashing' (LPH) we introduce and study in 
this Letter a new type of percolation model. Borrowing 
a name from the models of activated flow through traps 
Q we call this a 'drop-push' model since it has two basic 
features: a sequential 'dropping' of particles on a sub- 
strate followed by the transport of the dropped particles 
via a 'pushing' mechanism caused by the local hardcore 
repulsion between particles on the substrate. Unlike in 
the ordinary percolation, we show that the dynamics of 
the drop-push model generates nontrivial spatial correla- 
tions between sites. Our exact solution in one dimension 
shows that the critical exponents associated with the per- 
colation transition in the drop-push model are different 
from those of the ordinary percolation. As an additional 
bonus, our approach also rederives, in a straightforward 
way, a key result on the cost function in the LPH al- 
gorithm obtained recently by computer scientists using 
more involved combinatorial techniques . Our model 
is also easily generalizable to higher dimensions. 

The LPH algorithm was originally introduced by 
Knuth H and has remained popular in computer science 
due to its simplicity, efficiency and general applicability 
M. It can be described as follows: Consider M items X\, 
X2, • ■ ■ , xm to be placed sequentially into a table with L 
cells labelled 1, 2, . . ., L. For each item a;,-, a cell with 
label hi € {1,2, ... ,L} is selected. The label hi is called 
the hash address and is usually chosen randomly from 
the set {1,2, ... ,L}. The item Xi is inserted at the hi-th 
cell provided it is empty. Otherwise one tries cells hi + 1, 
hi + 2, etc. until an empty cell is found (the locations 
of the cells are interpreted modulo L) where the item X; 
is finally inserted. From the computer science point of 
view, the object of interest is the cost function S(M, L) 
defined as the total number of unsuccessful probes en- 
countered in inserting the M items into a table of size L 
[^],[|. This cost function measures the construction cost 
of the table as well as the time to search for an item later 



H||. The statistics of this cost function were found to 
be very different in the sparse table (when the density 
M /L « 1) compared to the full table (when M/L — > 1) 
0. We show below that, when interpreted as an inter- 
acting particle system, this crossover from the sparse to 
full table corresponds precisely to a percolation transi- 
tion which belongs to a different universality class than 
that of the ordinary site or bond percolation. 
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FIG. 1. The dropping and the subsequent hopping moves 
in the drop-push model. Also shown a particle cluster of size 
7 and a hole cluster of size 5. 



In our equivalent drop-push model (Fig. 1), we inter- 
pret the table as a lattice of size L with periodic bound- 
ary conditions. The cells are the lattice sites and each 
site can contain at most one particle. One starts with 
an empty lattice. At each step a particle is dropped into 
a randomly selected site. If the target site is empty, the 
incoming particle occupies it. If however the site was al- 
ready occupied, the particle keeps hopping to the right 
until it finds an empty site where it then stays (see Fig. 
1). One then repeats the same procedure with the next 
particle and so on. The dropping process in this model 
is similar to that of random sequential adsorption (RSA) 
and the car parking processes M, however the adsorp- 
tion mechanism in the drop-push model is quite different 
from the usual RSA models. In the car parking language, 
in this drop-push model, a new car arrives at a random 
spot on a one way lane and moves forward till it finds an 
empty parking spot. Unlike usual car parking models, 
here a car always manages to find a spot and thus the 
system never gets stuck in a jammed state. Note that al- 
though we have defined the hopping of the particle to be 
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unidirectional, one can also consider an unbiased version 
where a particle, if dropped onto an occupied cluster, 
performs an unbiased random walk on the occupied clus- 
ter till it finds an empty site where it is then attached. 
We show below that in one dimension, the results are 
independent of this bias. 

Each addition of a new particle corresponds to incre- 
menting the density t = M/L by the amount At = l/L. 
Thus in the thermodynamic limit L — ► oo, the den- 
sity t becomes a continuous 'time' like variable that in- 
creases monotonically from t = (empty lattice) to t = 1 
(full lattice). For convenience, we will henceforth refer 
to the density t as 'time' with < t < 1. The cen- 
tral objects of our analysis are P n (t) and Q n {t) denot- 
ing respectively the number (per unit length) of parti- 
cle and hole clusters of size n at time t (see Fig. 1). 
The total number of particle (hole) clusters is denoted 
by N(t) = £„P„(t) = £„<?»(*)■ Note that as t in- 
creases continuously from to 1, one expects that N(i) 
starting from N(0) = should increase initially but will 
eventually decrease to again as t — » 1 when the lattice 
is nearly full. Thus N(t) has an interesting nonmono- 
tonic behavior in < t < 1 with a maximum (when the 
system has the largest number of clusters) at an inter- 
mediate time t* (see the inset of Fig. 2). We also note 
the sum rules corresponding to the particle and the hole 
densities: £ nP n (t) = t and nQ n (t) = (1 - t). 




FIG. 2. The exact solutions for the normalized cluster den- 
sities p„ — Pn/N (filled triangles) and q n — Q„/N (filled 
squares) are compared with numerical results (circles and di- 
amonds respectively) obtained via the Monte Carlo simula- 
tion of the drop-push model on a lattice with L = 100000 
at the filling factor t = 0.5. The inset shows a plot of 
the numerical domain density N(t) and the exact result 
N(t) = (1 — t)(l — e _t ). The two curves are indistinguish- 
able. 



The key observation that allows the exact solution in 
one dimension is the fact that the 'cluster' or 'domain' 
densities P n 's and Q„'s are statistically independent at 
all times < t < 1. This follows from the fact that 
the adsorption of a new particle at the edge of a particle 
(hole) cluster does not introduce correlations between the 
adjacent particle and hole clusters. A rigorous justifica- 
tion of this fact, details of which will be published else- 
where ||, follows from the observation that the domain 
walls in the drop-push model can be viewed as the zero 
crossings of a Markov process in space (at a fixed time t). 
Thus if one labels a configuration C by the set {rii} where 
rij's denote the lengths of the alternate particle and hole 
clusters, then the probability of C is given by the prod- 
uct measure, Prob[C,i] oc P ni (t)Q n2 (t)P n3 (t)Q n4 (t) . . .. 
In other words, the independent interval approximation 
(IIA) is exact in this model. The next step is to write 
down the exact rate equation of evolution of P n 's and 
QnS by accounting for all the gain and loss terms due to 
the addition of a new particle and exploiting the factor- 
ization property of the probabilities. Similar types of IIA 
equations have been used in one dimensional coarsening 
problems jjj. The rate equation for the Q„'s turn out to 
be simple 
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valid for all n > 1. The first term denotes the loss of a 
hole cluster of size n due to an adsorption of a particle 
by direct hit at any of the available n sites and also due 
to a hit at any occupied site of the neighboring particle 
cluster to the left which then transports the particle to 
the leftmost site of the hole cluster. The latter happens 
with rate J2k kPk/N = t/N. The second term denotes a 
gain in Q n due to a direct hit inside a hole cluster of size 
bigger than n. The factor 2 denotes that there are only 
two places available for the incoming particle in order to 
generate a hole cluster of size n from a bigger hole cluster. 
The third term indicates a gain of a hole cluster of size n 
from that of size n+l due to a particle adsorption at its 
left edge. One can similarly write down the equations for 
the Pn's though they turn out to be trickier. Omitting 
the details M , we present only the final results 
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Although the above equations are written down for the 
unidirectional version of the model, a careful analysis 
shows that they remain unchanged for the general case 
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in which the dropped particle moves to the right with 
probability p or to the left with probability (1 — p) S. 

As an important consistency check, one can verify that 
both the Eqs. (|l|) and (^) satisfy the respective sum rules 

n Qn = (1 — t) and J2 n = O ne Can a l S0 wr h c 

down the evolution equation for the total domain density 
via direct inspection of the process, 

f = -2N-± Ql + l- t . (3) 

It is easy to check that both Eqs. (|l|) and (||) when 
summed over n satisfy Eq. (|3|) thus providing yet an- 
other useful consistency check. 

We note that Eq. (Q) for the Q n 's does not involve 
the P„'s, but is however implicitly nonlinear due to the 
occurrence of TV = J\ Q n on the right hand side. How- 
ever Eq. (Q) admits a very simple pure exponential so- 
lution as found in many RSA models M. The ansatz 
Q n (t) = A(t) exp[-nB(t)] satisfies Eq. (|) for all n > 1 
with the choice A(t) = 2(1 - i)[cosh(t) - 1] and B(i) = t. 
One then gets the exact solution for the domain den- 
sity N(t) = X)„ Qn(t) = (l-t)(l- e"*) which is non- 
monotonic in < t < 1 and is asymmetric about its 
unique maximum at t* = 0.4428 . . . (see Fig. 2). Note 
that in the ordinary site percolation in 1-d with occupa- 
tion probability t, the domain density is simply given by 
N(t) = t(l — t) which is symmetric about the maximum 
at t* = 1/2. 

We next substitute the exact solution for the Q„'s into 
Eq. (||) and first solve for n = 1 and n = 2. The exact so- 
lutions Pi{t) = t(l-t)e~ 2t and P 2 (t) = 3t 2 (l - t)e~ 3t /2 
suggest the ansatz: P n (t) — a n (l — t)t n e~( n+lS>t . Indeed 
this ansatz solves Eq. (|J) provided the o n 's satisfy the 
nonlinear recursion relation, 

n-2 

na n = (n + l)a n -i + ^J(« + l)oiO n -i-i, (4) 

i=l 

starting with a\ = 1. The first few values are a 2 = 3/2, 
a 3 = 8/3, a 4 = 125/24 etc. From Eq. it fol- 

lows that the generating function a(z) — _■, a n z n 
is given by a(z) = — 1 + T(z)/z where T{z) is the 
well known tree function given by the solution of the 
equation T(z) exp[— T(z)] — z. It is easy to see that 
the function T(z) has a singularity at z = 1/e where 
T(z) l — y/2(l — ez)+. . .. This tree function appears in 
problems related to the counting of rooted labelled trees 
fiof with various applications in computer science ||]^] 
as well as many physics problems such as aggregation 
models Jll|] and the classical hard sphere fluid in large 
dimensions |^2[ . Using the known properties of T(z), we 
then get a n = (n + l)" _1 /?i! for all n > 1. We have also 
checked that the exact solutions 

Q„(t) = (l-t)(e t -l) 2 e-("+ 1 ) t (5) 

P n (t) = (1 - t)t n e -in + l)t {n + lY- 1 



match perfectly with the numerical results obtained via 
the Monte Carlo simulation of the drop-push model (see 
Fig. 2). 

Clearly at t = 1, there is only one infinite parti- 
cle cluster and the system percolates. We now an- 
alyze the scaling behavior of the cluster distributions 
near the critical point t = 1. From Eq. (^), we find 
that for large n, P n {t) w ^=^- n~ e exp[— n/n*{t)\ where 
n*(t) = l/(t — 1 — log t) and the Parisi-Sourlas exponent 
]13[ is given exactly by 6 = 3/2. Note that for the or- 
dinary percolation in one dimension, P n (t) = (1 — t) 2 t n 
indicating 9 — 0. In the limit t — > 1, the typical clus- 
ter size diverges as n*(t) « 2(1 — t)~ 2 and one ob- 
tains the Stauffer scaling form M for the cluster size, 
P n (t) ~ n~ T f[n(\ — t) a ] with the Fisher exponent r = 2 
(as in the ordinary percolation) and a — 2 (in contrast 
to the ordinary percolation where a = 1). The exact 
scaling function here f(z) = \/^"e _z / 2 also differs from 
that of the ordinary percolation where f(z) — z 2 e~ z [EJ. 
Consequently the susceptibility exponent ffl given by the 
scaling relation 7 = (3 — t)/<j also differs in the two mod- 
els. For the drop-push model 7 = 1/2, where as 7 = 1 
for the ordinary percolation. 

We now turn to the correlation function G n (t) denot- 
ing the probability that two occupied sites separated by 
a distance n belong to the same cluster at time t. For 
convenience we introduce the binary variable cr^ such 
that <7{ = 1 if the site i is occupied and <Ji = oth- 
erwise. Then by definition G n = (<Ji<j i+ i . . . o i+n ). We 
also note that by definition the particle cluster density 
P n = (aicr l+1 . . . <j i+n a l+n+1 ) where a* = 1 - a^. As 
a consequence one obtains the exact relation, G n +i + 
G„_i — 2G n — P n . Using the exact scaling form of P n 
and integrating twice with respect to n we find that in the 
scaling limit t — ► 1, n — > 00, G n (t) w g[n(l — t)"] where 
the correlation length exponent v — 2 and the exact scal- 
ing function g(z) = (1 + z)crfc( v /|) - y^-e -2 / 2 . These 

results should be compared to those for the ordinary per- 
colation with occupation density t where G n (t) = t n in- 
dicating that v — 1 and g(z) = e~ z trivially. 

To elucidate the nontrivial spatial correlations in the 
drop-push model we have also computed the conven- 
tional two point correlation function, C n (t) = ((Tj(Tj-|_ n ) — 
(<Ti) ((Ti+n) , the connected part of the joint probability 
that two sites at distance n are both occupied. Note that 
for the ordinary percolation, C n (t) — for all n > 1 and 
< t < 1. In contrast, we show that C n [t) is nontrivial 
in the drop-push model. In order to compute it, we add 
up the possibilities that there may be no holes between 
the two sites, or maybe only one hole cluster, or two hole 
clusters etc. This method of expressing the correlation 
function in terms of the interval size distributions was 
used before in other contexts |l4|,[l5|]. Omitting details 
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we present only the final expression for the generating 
function C{z) = £ n C n (t)z n 



C{z) 



1 - 1 

1-2 



Z*- 1 



t(l-z) 



(J) 



where T(z) is the tree function defined earlier. Using 
the properties of the tree function it is straightforward, 
though somewhat tedious ||, to derive the asymptotic 
scaling properties of C n (t). We find after some algebra 
that in the scaling limit t — ► 1, n — ► oc but keeping 
z = n(l-t) 2 fixed, C n (t) fa (1 -i) 2 F[n(l -t) 2 l where the 
exact scaling function F(z) = ,} e~ 

« 1/ " 



The function 

1 „-z/2 




-z/2 



ierfc( v ^ 
and F(2 



as z 



v^z 3 / 2 

We now turn to the total cost function S(t) in the 
LPH algorithm, i.e. in the drop-push model with uni- 
directional transport. The cost AS(t) to insert a new 
particle at time t is precisely the expected number of 
unsuccessful probes, i.e., the expected number of hops 
that the new particle undergoes before getting adsorbed 
into the right edge of a particle cluster. Consider a 
particle cluster of size k at time t. The incoming par- 
ticle can drop anywhere on this cluster and the num- 
ber of hops is simply the distance of the target site 
from the right edge of the cluster. Noting that the 
cluster size can vary from 1 to infinity, one then gets 

*s(t) = YZi (Ej=ii)*M*) = hEZiKk + mit). 

Using the exact result for Pk(t) from Eq. (g) we get 
AS(t) = t{2 - t)/[2(l - t) 2 ]. The total cost function 
is then given by S(t) = f Q AS(t')dt' = t 2 /[2(l - t)], in 
agreement with the result derived by the computer scien- 
tists using rather involved combinatorial techniques plH . 
More details on the statistics of this cost function for fi- 
nite size tables have also been derived recently (|U • 

Finally this drop-push model can be easily generalized 
to higher dimensions. Let us consider, for simplicity, the 
the unbiased version. One starts with an empty lattice of 
linear size L in d-dimensions with periodic boundary con- 
ditions. At each step one drops a particle at a randomly 
chosen site. The incoming particle occupies the target 
site provided it was empty. Otherwise the dropped par- 
ticle performs a random walk starting at the target site 
until it finds an empty site and sticks there. One then 
adds another particle and the process continues. It fol- 
lows from the simple electrostatic analogy that when a 
particle drops onto an occupied cluster, it has equal prob- 
ability to subsequently stick to any of the surface sites of 
this cluster. In this sense this model to similar to the 
celebrated Eden model However, unlike the Eden 

model, here one can have many different seed sites from 
which a new cluster can grow. Also the probability that a 
given cluster will grow by absorbing a new particle is pro- 
portional to the volume of the cluster in the drop-push 
model. Note that this model is also different from the pre- 
viously studied cooperative adsorption models |17|, the 



cluster-cluster aggregation models |16[ and the random 
dynamical percolation model |18|. It is clear that there 
will be a critical density t = t c < 1 at which the particle 
clusters start to percolate. An outstanding question is 
whether this percolation transition in higher dimensions, 
as in the 1-d case, belongs to a different universality class 
than that of the ordinary site percolation. Our prelim- 
inary numerical simulations in 2-d indicate that indeed 
this may be the case ||] . We defer these results and other 
details for a future communication. 

We thank M. Barma and D. Dhar for useful discus- 
sions. 
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